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We present the results of a study of the charmless vector- vector decay B + —* p + K*°, based on 
253 fb _1 of data collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. 
We obtain the branching fraction B{B+ -> p + K*°) = (8.9 ± 1.7(stat) ± 1.2(syst)) x 10" 6 . We also 
perform a helicity analysis of the p and K* vector mesons, and obtain the longitudinal polarization 
fraction f L {B+ -> p+ K*°) = 0.43 ± 0.11(stat)t°;^ (syst). 

PACS numbers: 13.25.Hw, 14.40.Nd. 



Naive factorization in the Standard Model (SM) pre- 
dicts that the longitudinal polarization fraction {Jl) hi 
B meson decays to light vector-vector (VV) hnal states 
is close to unity 0. In the tree dominated B + — ► p + p a 
and B° — > p + p~ decays, this prediction has been con- 
firmed In contrast, for the pure b — > s penguin 
B — > 4>K* decay, Belle and BaBar Q have found 
that the longitudinal and transverse polarization frac- 
tions are comparable, which is in disagreement with the 
factorization expectation. Possible explanations for this 
discrepancy include enhanced non-factorizable contribu- 
tions such as penguin annihilation [3, large 577(3) break- 
ing in form factors or new physics 7, 8J. It is therefore 
important to perform polarization measurements in other 
W modes, in particular, in the pure penguin b — ► sdd 
decay B+ -> p+K*°. 

In this paper, we present the results of a study of 
B+ -> p+K*° decays with a 253 fb" 1 data sample 
containing 275 x 10 6 B meson pairs collected with 
the Belle detector at the KEKB asymmetric-energy 
e + e~ collider [T^j operating at the T(4S) resonance 
(Vs = 10.58 GeV). The production rates for B+B- and 
B°B° pairs are assumed to be equal. 

The Belle detector is a large solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside of the coil is instru- 
mented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [llj . 



We select B + — > p + K*° candidate events by combining 
three charged tracks (two oppositely charged pions and 
one kaon) and one neutral pion. Each charged track is 
required to have a transverse momentum px > 0.1 GeV /c 
and to have an origin within 0.2 cm in the radial direction 
and 5 cm along the beam direction of the interaction 
point (IP). 

Particle identification likelihoods for the pion and kaon 
hypotheses are calculated by combining information from 
the TOF and ACC systems with dE/dx measurements 
in the CDC. To identify kaons, we require the kaon like- 
lihood ratio, Lk/{Lk + L v ), to be greater than 0.6. To 
identify pions, we require Lk/{Lr + L^) to be less than 
0.4. The efficiency for this selection is 86% for kaons and 
89% for pions, with corresponding tt/K misidentification 
rates of 8% and 10%. In addition, charged tracks are re- 
jected if they are consistent with an electron hypothesis. 

Candidate 7r° mesons are reconstructed from pairs of 
photons that have an invariant mass in the range 0.1178 
- 0.1502 GeV/c 2 , corresponding to a window of ±3er 
around the nominal ir° mass. The photons are assumed 
to originate from the IP. The energy of each photon in the 
laboratory frame is required to be greater than 50 MeV 
for the ECL barrel region (32° < 9 < 129°) and 100 
MeV for the ECL endcap regions (17° < 9 < 32° and 
129° < 6 < 150°), where 9 denotes the polar angle of the 
photon with respect to the beam line. The ir° candidates 
are kinematically constrained to the nominal ir° mass. In 
order to reduce the combinatorial background, we only 
accept 7T° candidates with momenta p^o > 0.40 GeV/c 
in the e + e~ center-of-mass (CM) system. 

Candidate p + mesons are reconstructed via their 
p + — * 7T + 7r° decay, and the 7r + 7r° pairs are required to 
have an invariant mass in the region 0.62 GeV/c 2 < 
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M (tt+tt ) < 0.92 GeV/c 2 . Candidate K*° mesons are 
selected from the K*° — > if + 7r~ decay with an invariant 
mass 0.83 GeV/c 2 < M(K + n~) < 0.97 GeV/c 2 . 

To isolate the signal, we form the beam-constrained 
mass Mbc = yj £^ cam — , and the energy difference 
AE = Eb — -Ebcam, where i?beam is the CM beam en- 
ergy, and pb and Eb are the CM momentum and energy, 
respectively, of the B candidate. The AE distribution 
has a tail on the lower side caused by incomplete lon- 
gitudinal containment of electromagnetic showers in the 
CsI(Tl) crystals. We accept events in the region M\, c > 
5.2 GeV/c 2 and -0.3 GcV < AE < 0.3 GcV, and define 
a signal region in M^ c and AE as 5.27 GeV/c 2 < M^ c < 
5.29 GeV/c 2 and -0.10 GeV < AE < 0.06 GeV respec- 
tively. These requirements correspond to approximately 
±3 a for both quantities. 

The continuum process e + e~ — > qq (q = u,d,s,c) is 
the main source of background and must be strongly sup- 
pressed. One method of discriminating the signal from 
continuum is based on the event topology, which tends 
to be isotropic for BB events and jet-like for qq events. 
Another discriminating characteristic is 8b, the CM po- 
lar angle of the B flight direction. B mesons are pro- 
duced with a 1 — cos 2 9b distribution while continuum 
background events tend to be uniform in cos 6b ■ The 
displacement along the beam direction between the sig- 
nal B vertex and that of the other B, Az, also provides 
separation. For B events, the average value of Az is 
approximately 200 /im, while continuum events have a 
common vertex. Additional discrimination is provided 
by the 6-flavor tagging algorithm 0] developed for time- 
dependent analysis at Belle. The flavor tagging proce- 
dure yields two outputs: q (= ±1), which indicate the 
flavor of the tagging B, and r, which ranges from to 
1 , is a measure of the likelihood that the b flavor of the 
accompanying B meson is correctly assigned. For signal 
events, q is more likely consistent with the opposite of the 
charge of signal B; there is no correlation for continuum 
events. Events with high values of r are well-tagged and 
are less likely to originate from continuum production. 
Thus, the quantity q ■ r ■ C'b, where C'b is the charge of 
the signal B, can be used to discriminate against contin- 
uum events. 

We use Monte Carlo (MC) simulated signal and data 
sideband (defined as 5.2 GeV/c 2 < M bc < 5.26 GeV/c 2 ) 
events to form a Fisher discriminant based on a set of 
modified Fox- Wolfram moments 0] that are confirmed 
to be uncorrelated with Mb c , AE and variables consid- 
ered later in the analysis. Probability density functions 
(PDFs) derived from the Fisher discriminant, the cos 9b 
distributions and the Az distributions are multiplied to 
form likelihood functions for signal (jC s ) and continuum 
background {Cqq)', these are combined into a likelihood 
ratio TZ S = £ S /(C S + C q q). We achieve background sup- 
pression by imposing g-r-C^-dependent 1Z S requirements, 
which are determined by optimizing the figure of merit, 



SyV S + B, where S (B) is the number of signal (back- 
ground) events in the signal region. A branching fraction 
of B(B + — > p + K *°) = 1 x 10~ 5 is assumed. This require- 
ment removes 99.3% of the continuum background while 
retaining 41% of the B+ -> p+K*° events. The MC- 
determined efficiency with all selection criteria imposed 
is 2.7% for longitudinal polarization (Aq) and 4.0% for 
transverse polarization (A±). 

The fraction of multiple candidates in the signal region 
for signal MC is 3.6% for the Aq helicity state and 1.7% 
for the A± state. We allow multiple candidates in this 
analysis. 

To investigate backgrounds from b — > c decays, we use 
a sample of BB MC events corresponding to an inte- 
grated luminosity of 412 fb _1 . We find a contribution 
from B + — > U° '(K + tt~ tt )tt + decays in the p or K* 
sideband region and require \M(Kirir°) — M D o\ > 0.050 
GeV/c 2 to veto these events. This requirement does not 
remove any B + — ► p + K*° events. Among the charmless 
B decays, potential backgrounds arise from B + — > a\K + , 
B+ -> p+K^(U30), non-resonant B+ -> p+K+n- and 
B + — ► K*°tt + tt°. We separate signal from these back- 
grounds by fitting the p and K* invariant mass distribu- 
tions. 

We extract the signal yield by applying an extended 
unbinned maximum-likelihood fit to the two-dimensional 
Mbc-AE distribution. The fit includes components for 
signal plus backgrounds from continuum events and b — > 
c decays. The PDFs for signal and b — > c decay are 
modeled by smoothed two-dimensional histograms ob- 
tained from large MC samples. The signal PDF is ad- 
justed to account for small differences observed between 
data and MC for a high-statistics mode containing ir° 
mesons, B + — > D°(K + ir~ir a )ir + . The continuum PDF 
is described by a product of a threshold (ARGUS) func- 
tion 0] for Mbc and a first-order polynomial for AE, 
with shape parameters allowed to vary. All normaliza- 
tions are allowed to float. Figure ^ shows the final event 
sample and the fit results. The five-parameter (three nor- 
malizations plus two shape parameters for continuum) fit 
yields 134.8 ± 16.9 B+ -> K+ir-ir+ir° events. 

We further distinguish the p + K*° signal from non- 
resonant decays such as B + — > p + K + n~ or B + — > 
^* 7r+7r° by fitting the M (tt+tt ) and M(K+ir~) invari- 
ant mass distributions. The signal yields obtained from 
the M bc - AS fit for different M(tttt) and M(Kn) bins 
are plotted in Fig. where the M(7ttt) distribution is 
for events in the K* region (0.83 GeV/c 2 < M(Kw) < 
0.97 GeV/c 2 ) and the M(Kir) distribution is for events 
in the p region (0.62 GeV/c 2 < M(tttt) < 0.92 GeV/c 2 ). 
We perform separate x 2 fits to the M(tt7t) or M(Kir) 
distributions. Each fit includes components for signal 
and non-resonant background. The signal p and K* 
PDFs arc modeled by relativistic P-wave Breit-Wigner 
functions with means and widths fixed at their known 
values 0|; the PDFs are convolved with a Gaussian of 
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FIG. 1: Projections of M^ c for events in the AE signal region 
(left), and projection of AE in the Mt, c signal region (right). 
The solid curves show the results of the fit. The hatched 
histograms represent the continuum background. The sum of 
the b — > c and continuum background component is shown as 
dot-dashed lines. 



fj = 5.3 MeV, which is obtained by fitting the D a (K~ir + ) 
invariant mass, to account for the detector resolution. 
The non-resonant component is represented by a thresh- 
old function with parameters determined from MC events 
where the final states are distributed uniformly over 
phase space. The M(irir) mass fit gives 125.4 ± 15.8 p 
and —0.3 ± 3.0 non-resonant K*°Tr + ir° events in the p 
mass region. In the M(Kir) fit, we find 85.4 ± 16.1 pK* 
signal and 28.8 ±4.1 non-resonant events in the K* mass 
region. The statistical significance of the signal, defined 
as y/Xo ~ Xmin> where Xmin is the X 2 value at the best-fit 
signal yield and Xo is the value with the K* a signal yield 
set to zero, is 5.3 a (5.2 a with the inclusion of system- 
atics). The contribution from non-resonant p + K + n~ is 
significant and is taken into account in both the branch- 
ing fraction and polarization determinations, while we 
neglect the non-resonant K*°tt + tt° contribution. 
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FIG. 2: Signal yields obtained from the M^ C -AE distribution 
in bins of M(tv + tt°) (left) for events in the K*° region and 
in bins of M(K + tt~) (right) for events in the p region. Solid 
curves show the results of the fit. Hatched histograms are for 
the non-resonant component. 

We use the p + — > 7r + 7r° and K*° — » K + tt~ helicity- 
angle (9 P , 9k*) distributions to determine the relative 



strengths of \A \ 2 and |^4±| 2 - Here 9 P (9k*) is the angle 
between an axis anti-parallel to the B flight direction and 
the 7t + (K + ) flight direction in the p (K*) rest frame. For 
the longitudinal polarization case, the distribution is pro- 
portional to cos 2 9 p cos 2 9k* , and for the transverse po- 
larization case, it is proportional to sin 2 9 P sin 2 9k * |16| . 
Figure [3] shows the signal yields obtained from M\> C -/S.E 
fits in bins of the cosine of the helicity angle for p and K* . 
We perform a binned simultaneous \ 2 fit to the p and K* 
helicity-angle distributions. The fit includes components 
for signal and non-resonant pKn. PDFs for signal Aq and 
A± helicity states are determined from the MC simula- 
tion. The helicity- angle distribution for data in the high 
M (Kir) sideband region 1.1 GcV/c 2 < M(Kn) < 1.5 
GeV/c 2 , where pKir events dominate, is consistent with 
a cos 2 9-\ike cos 9 P and a flat cos 9k* distribution. Thus, 
we assume an S'-wave Kir system and model the non- 
resonant B — > pKir PDF based on the MC simulation. 
The fraction of the non-resonant component is fixed at 
the values obtained from the K* mass fit. The two- 
parameter (normalizations for Aq and A±) fit result de- 
viates from 100% longitudinal polarization with a signif- 
icance of 4.9ct (4.3cr including systematic uncertainties). 
The significance is defined as \Jxq~ Xmhv wnere Xmin is 
the x 2 value at the best-fit and Xo is the value with the 
longitudinal polarization fraction set to 100%. 

The largest uncertainties in the polarization measure- 
ment are due to uncertainties in the non-resonant pKir 
PDF, potential scalar-pseudoscalar (S-P) interference, 
and the non-resonant fraction. We assign a 1q°' 3 % sys- 
tematic error for the non-resonant PDF. This uncer- 
tainty is estimated by adding a 1/3 flat component to 
the p helicity PDF for non-resonant pKir in the helic- 
ity fit. Interference of the longitudinal amplitude Aq 
with the S'-wave (Kir) system introduces a term with a 
2e % ^\A P K-n\ cos 9k* dependence, where A<j> is the phase 
difference and \A p k-k\ is the amplitude of the B — > pKir 
decay. The S-P wave interference disappears in the cos 9 P 
distribution, which is integrated over cos9k*] however it 
remains in the cos 9k* distribution. We include an ad- 
ditional linear function for the interference term in the 
cos9k* helicity, and redo the x 2 fit- The resulting small 
change in /l, 0.5%, is assigned as the systematic un- 
certainty for the S-P interference. A ^4'°% systematic 
error is assigned for the uncertainty in the fraction of 
non-resonant pKir, obtained by varying the non-resonant 
fraction by ±1(7. Adding the various systematic error 
contributions in quadrature, we obtain the longitudinal 
polarization fraction in B + — ► p + K*° decays 

f L (B+ -> p+K*°) = 0.43±0.11(stat)±8:8|(syst). 

To calculate the B + — > p + K*° branching fraction, 
we use the M(Ktt) invariant mass fit result and MC- 
determined efficiencies weighted by the measured polar- 
ization components. We consider systematic errors in the 



5 



60 
50 
40 

© 

g 30 
? 20 
10 




60 r- 

50 ~ r 

•* 40 - 
© 

? 20 ^ 
10 - 



-1 



-0.5 

cos8„ 



o 
-l 



-0.5 

cos8„ 



FIG. 3: Fit to background-subtracted helicity distributions. 
The solid histograms show the results. The dot-dashed 
(dashed) histograms are the Ao (A±) component of the fit; 
the dotted histograms are for non-resonant pK-n. The low 
event yield near cos# p = 1 is due to the p n o > 0.4GeV/c 
requirement. 



branching fraction that are caused by uncertainties in 
the efficiencies of track finding, particle identification, 7r° 
reconstruction, continuum suppression, fitting, polariza- 
tion fraction. We assign an error of 1.1% per track for the 
uncertainty in the tracking efficiency. This uncertainty 
is obtained from a study of partially reconstructed D* 
decays. We also assign an uncertainty of 0.7% per track 
on the particle identification efficiency, based on a study 



of kinematically selected D* + — > D 



D° 



K-TT+ 



decay. A 4.0% systematic error for the uncertainty in 
the 7T° detection efficiency is determined from data-MC 
comparisons of r\ — ► 7r°7r°7r with rj — > ir + ir~ir and 
V ~ * 77- A 4.5% systematic error for continuum suppres- 
sion is estimated from studying the process B + — > D°tt + , 
~D° -> K+tt-tt . A -2.0%/ + 1.7% systematic error as- 
sociated with fits is obtained by shifting each parameter 
by ±1 a. A 6.7% systematic error for the uncertainty in 
the b — ► c background PDF is obtained by changing the 
PDF parameterization. A — 4.2%/+4.4% error due to the 
uncertainty in the fraction of longitudinal polarization is 
obtained by varying fa by its errors. The uncertainty in 
non-resonant K*irir background gives a contribution of 
-2.2%/ + 0% in addition to -3.0%/ + 2.3% error from 
uncertainties in the background from other rare B de- 
cays. A 7.1% error for possible bias in the x 2 n t 113 is 
obtained from a MC study. A 1.1% error for the uncer- 
tainty in the number of BB events in the data sample 
is also included. The quadratic sum of all of these er- 
rors is taken as the total systematic error. We obtain the 
branching fraction 



B(B^ 



f iT°) = 



3.9 ± 1.7(stat) ± 1.2(syst)) x 10" 



We measure the branching fraction to be 
(8.9 ± 1.7(stat) ± 1.2(syst)) x 10" 6 . We also perform 
a helicity analysis and find a substantial transversely 
polarized fraction with a statistical significance of 4.9 a. 
The longitudinal polarization fraction fa measured is 
similar to the surprisingly low value found in b — > sss 
decays B — ► <ftK* . 
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In summary, we have observed the B + — ► p + K* 
decay with a statistical significance of 5.3 a. 
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